Quantum conditional operator and a criterion for separability 
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We analyze the properties ol the conditional ampli tude operator, the quantum analog of the 



conditional probability which has been introduced in [}uant-ph/9512022|. The spectrum of the 
conditional operator characterizing a quantum bipartite system is invariant under local unitary 
transformations and reflects its inseparability. More specifically, it is shown that the conditional 
amplitude operator of a separable state cannot have an eigenvalue exceeding 1, which results in a 
necessary condition for separability. This leads us to consider a related separability criterion based 
on the positive map F : p — » (Trp) — p, where p is an Hermitian operator. Any separable state is 
£ — , mapped by the tensor product of this map and the identity into a non-negative operator, which 

^\ • provides a simple necessary condition for separability. In the special case where one subsystem is a 

ON \ quantum bit, F reduces to time-reversal, so that this separability condition is equivalent to partial 

transposition. It is therefore also sufficient for 2x2 and 2x3 systems. Finally, a simple connection 
between this map and complex conjugation in the "magic" basis is displayed. 

O ! PACS numbers: 03.65.Bz, 89.70.+C KRL preprint MAP-217 



I. INTRODUCTION 



The state of a quantum bipartite system AB is described as separable (or classically correlated) if it can be obtained 
by having both parties A and B preparing their subsystem according to some common instructions (see, e.g., Jl|,||). 
Mathematically, this means that the density operator p characterizing the state of the bipartite system can be written 
£~^. . as a convex sum of product states, that is 



o 
o 



where the weights Wi satisfy Wj = 1 and < Wi < 1. The w^s can be viewed as the probability distribution of 
—j ■ a classical random variable that is known to both parties A and B and used by them to prepare their subsystem. 

i (A} (£3) 

Namely, if the subsystem A (and B) is prepared in state p\ ' (and p\ ') when the classical variable takes on value i, 
the state of the joint system AB is given by Eq. ([!]). A separable state p satisfies several interesting properties. The 
, joint statistics of any pair of local observables Oa and Ob (measured separately on each subsystem) can be described 
classically, based on an underlying global "hidden" variable. For example, the quantum expectation value of the 
product OaOb is given by 



Tr[p(0 A ®0 B )} =J2 w i( a )i( b )i ( 2 ) 

i 

where (a)i = Tv[pl Oa] and (b)i = Ty[p\ B ^ Ob\- In other words, the joint statistics of Oa and Ob can be understood 
classically, by assuming that the local statistics of the outcomes can be described separately for each p^' and p\ , 
and that the correlations originate from a hidden variable i distributed according to u>i. Moreover, a separable system 
always satisfies Bell's inequalities (the converse is not true), so that the latter represent a necessary condition for 
separability (see, e.g., Q). Note that any joint probability distribution can be represented as a convex combination 
of product distributions, so that classical probabilities are always separable in the above sense. 

The decomposition of a separable state p into a convex mixture of product states is not unique in general, but the 
fact that p is separable implies that there must exist at least one such decomposition. If no such decomposition can be 
written, then p is termed inseparable, and it can be viewed as quantum correlated. Except for the special case where 
p describes a pure state, the distinction between separable and inseparable states appears to be an extraordinary 
difficult problem. More precisely, some mixed states can be "weakly" inseparable, in the sense that it is very hard 
to establish with certainty their inseparability. This is basically due to the difficulty of enumerating explicitly all 
the possible convex combinations of product states in order to detect that a state is actually inseparable. Still, it is 
possible to find some conditions that all separable states must satisfy, therefore allowing the detection of inseparability 
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when a state violates one such condition. The most common example of such a necessary condition for separability is 
the satisfaction of Bell's inequalities. A state that violates Bell's inequalities is inseparable, while a state satisfying 
them may be separable or weakly inseparable 0. 

More recently, a surprisingly simple necessary condition for separability has been discovered by Peres ^| , which has 
been shown by Horodecki et al. (^] to be strong enough to guarantee separability for bipartite systems of dimension 
2x2 and 2x3. If the state p is separable, then the operator obtained by applying a partial transposition with respect 
to subsystem A (or B) to p must be positive, that is 

p ta = (p t °Y>o (3) 

Thus, this criterion amounts to checking that all the eigenvalues of the partial transposition of p are non-negative, 
which must be so for all separable states. In Hilbert spaces of dimensions 2x2 and 2x3, this condition is actually 
sufficient, that is, it suffices for ruling out all inseparable states ||. In larger dimensions, however, it is provably 
not sufficient, in the sense that it does not detect some weakly inseparable states UK- A general necessary and 
sufficient condition for separability in arbitrary dimensions has been found by Horodecki et al. || , which states that 
p is separable if and only if the tensor product of any positive^] map (acting on ^4) and the identity (acting on B) 
maps p into a positive operator. Although very important in theory, this criterion is hardly more practical than the 
definition of separability itself since it involves the characterization of the set of all positive maps. It appears to be 
useful mainly for 2 x 2 and 2x3 bipartite systems, where such a general characterization has been found [||. 

In this paper, we focus on the connection between quantum non-separability and the conditional amplitude operator 
which has been introduced recently in the context of quantum information theory In Section II, we start 

by detailing the mathematical properties of such an operator (support, spectrum, connection with von Neumann 
entropies, etc.). We then derive a necessary condition for separability, based on the conditional von Neumann entropy 
and its underlying conditional amplitude operator. Namely, the eigenvalues of the latter operator cannot exceed 1 if 
the bipartite state is separable, as was conjectured in Refs. [^j-jf). Since the conditional von Neumann entropy can 
be negative only if the conditional amplitude operator admits an eigenvalue larger than 1, a related — and weaker — 
separability condition is that the conditional entropy is non- negative (see also Refs. and ||). This leads us to 

consider a positive map r : p — » (Trp) — p which gives rise to a simple necessary condition for separability in arbitrary 
dimensions. More specifically, it is shown in Section III that any separable state is mapped by the tensor product of 
r (acting on one subsystem, A) and the identity (acting on the other, B) into a non-negative operator. In the case 
where T is applied to a two-state system (quantum bit or spin-1/2 particle), this corresponds to the time-reversal 
operation applied on one system with respect to the other one. Since Peres' criterion has been shown to be unitarily 
equivalent to such a "local" time-reversal by Sanpera et al. ||, our separability criterion is simply equivalent to Peres' 
for 2 x n composite systems. Therefore, it also results in a sufficient condition for 2x2 and 2x3 systems, according to 
Ref. ||. It also has a very simple geometric representation in the Hilbert-Schmidt representation of the bipartite state. 
Finally, it appears that the map T is connected to the complex conjugation operation in the "magic" basis introduced 
by Hill and Wootters |]l0| . In Appendix A, we illustrate the separability condition based on T by applying it to several 
separable or inseparable states, and compare it to the separability condition based on partial transposition. 

II. CONDITIONAL "AMPLITUDE" OPERATOR 

Let us a consider a bipartite system AB, characterized by a density operator pab in the product Hilbert space 
H.AB = T~La <B> Ti-B- Each subsystem, A or B, is characterized by the reduced density operator pa — Ttb[pab] or 
Pb = Tta[pab], respectively. 



Definition 1: Define the conditional amplitude operator of A conditional on B as 

Pa\b = exp[logp AB - log(l A <8> Pb)} 

= limbY^CU®/^)- 1 /"]" (4) 

n^oo 

which is a positive semi-definite Hermitian operator in the joint Hilbert space Hab defined on the range of pab (see 
Lemma 1). 



A map is defined as positive if it maps positive operators into positive operators. 
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The second expression in Eq. (|j), based on the Trotter decomposition of pa\b, explicitly emphasizes that the 
conditional amplitude operator is the natural quantum analog of the conditional probability, p(a\b) — p(a, b) /p(b). As 
Pab and (1a ® Pb) -1 do not necessarily commute, the Trotter symmetrization is used here to obtain an Hermitian 
operator. 

Lemma 1: Ran(/?As) Q Ran(lA ® ps), where Ran(p) is the range of p.f] Thus, the support of pab is included in 
the support of \a® Pb, and, consequently, the conditional amplitude operator pa\b is well-defined on the support of 

PAB- 

We must prove that Ker(lA ® ps) Q Ker(pAs), where Ker(p) is the kernel of p, i.e., that any eigenvector of 
(1a <8> Pb) with zero eigenvalue is such that pabIV") — 0. First note that any such eigenvector can be written as a 
linear combination of states \<p) where 

\4>) = \a) ® \b) (5) 

and | a) is an arbitrary state vector in H.a while \b) is an eigenvector of ps with zero eigenvalue, i.e., pbW) = 0. Let 
us now consider the positive semi-definite operator p = (1a <S> Pb)pAB(XA <8> Pf>), with Pf, = \b)(b\. It is trivial to check 
that its partial trace over A vanishes, that is 

Tr a [p] = PbPBPb = (6) 

This results from the general relation 

Tt a [{1a <S> \b)pab] = X b Tt a [pab] (J) 

where Xg and pab are arbitrary operators in H.b and TLab respectively. Since p is positive semi-definite and traceless, 
we have p = 0. Thus, in particular, the expectation value of p in the state \<p) vanishes, 

(0|p|0)H0|pA B |0) = O (8) 

which in turn implies that Pab\<P) — since pab is positive semi-definite. As this is true for each term \<j>) in the 
superposition, we conclude that pab\^) = 0. □ 

Remark: Lemma 1 clearly implies that Ker(lA®ps)nRan(pAs) = 0- Thus, the subspace spanned by the eigenvectors 
with zero eigenvalue of (1a ® Pb) is disjoint from the support of pab, so that the definition of pa\b contains no 
singularities in the support of pab- Of course, there is a classical analog for probability distributions which ensures 
that p(a\b) = p(a,b)/p(b) is well defined if a, b are such that p(a,b) ^ 0. Indeed, if b is such that p(b) — 0, then 
p(a, b) — 0, Va. This is obvious since p(b) = J2 a p(a, b) and p(a, b) > 0. 

Definition 2: The conditional von Neumann entropy is defined using the joint density operator pab and the condi- 
tional amplitude operator pa\b as J|-|7j 

S(A\B) = -Tv[ P ABlog 2 pA lB } (9) 

in close analogy to the classical definition 

H(A\B) = -Y,P(a,b)log 2 p(a\b) (10) 

a,b 

Thus, S(A\B) corresponds to the quantum entropy of A conditional on B, and is mathematically well-defined as a 
consequence of Lemma 1. The trace in Eq. (^) can indeed be restricted to the support of pab, using the fact that the 
common eigenvectors \ip) with zero eigenvalue of both pab and (1a®Pb) yield a vanishing contribution to the entropy, 
as lim^^o log x) = 0. (The same argument is used in classical information theory to discard zero probabilities when 
calculating Shannon entropies.) 



2 For a linear and Hermitian operator p, the range is the subspace of the domain of p that is spanned by the eigenvectors 
corresponding to non-zero eigenvalues, i.e., it is the support of p. 
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Theorem 1: The definitions of pa\b and the conditional von Neumann entropy imply that S(A\B) — S(AB) — S(B), 
as for Shannon entropies. 

First, using Eqs. (||) and (||), we have 

S(A\B) = -Tr[ PAB log 2 PAB ] + Tr\p AB \og 2 {l A ® p B )\ (11) 

where the first term on the right-hand side is clearly equal to S(AB). In order to calculate the second term on the 
right-hand side of Eq. (|ll]), we write 

Tta[pab log 2 (lA ® Pb)] = TtaIpabO-A ® log 2 Pb)] 
= Tr aIpab] log 2 Pb 

= Pb log 2 Pb (12) 
where we have made use of Eq. (0). This implies that the second term on the right-hand side of Eq. (|li"| ) is 

Tv B [p B log 2 p B ] = -S(B) (13) 

resulting in S(A\B) = S(AB) - S(B). □ 

Lemma 2: The spectrum of the conditional amplitude operator pa\b is invariant under unitary transformations of 
the product form Ua®U b on pab- 

Let us consider the isomorphism 

Pas -> Pab = 

{U a ®U b )pab{U a ®uI) (14) 
We first calculate the partial trace of the joint density operator over A after this transformation, that is 

Pb = Tt a[p'ab] 

= Ty a [(Ua ® U b )pab{U\ ® Ul)] 

= Tr A [(l^ (g) U B ){U A ® 1b)pab{U\ ® 1 b )(1a ® C^)] 

= C/ s Tr A [((7A ® 1 b )pab(C^ ® lfl)]^ 

= C/bPbC/]? (15) 

where we have used Eq. (Q) and the basis invariance of the trace. This implies that the conditional amplitude operator 
transforms as 

Pa\b -> Pa\b = 

(Ua®U b ) P a\b{u\®U b ) (16) 

so that its spectrum is conserved under Ua ®!/b on Pab- Note that the classical analog of an U a ® Ub isomorphism 
corresponds to permuting the rows and columns of the joint probability distribution p(a, 6), so that the classical 
counterpart of Eq. (|l^) is straightforward. □ 

Remark: This Lemma suggests that the spectrum of Pa\b could be related to the separability of the state pab, since 
separability (or inseparability) is conserved under a U A ®Ub isomorphism. This will be examined later on. 

Corollary: The conditional von Neumann entropy S{A\B) is invariant under a unitary transformation of the product 
form Ua®Ub- 

This property results from the definition of S(A\ B), Eq. (^), together with Eq. (plf), or can be checked trivially 
from Theorem 1. 

Theorem 2: The operator oab = — ^ogPA\B = k>g(l A <S> Pb) — log pab is positive semi-definite if the quantum 
bipartite system characterized by pab is separable. 

Let us consider a separable bipartite system pab, i-e., a convex combination of product states: 

Pab = ^2 Wi (p^a ® w i tn w i — ^ an d < W{ < 1 (17) 
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where p^ and are states in TLa and TLb, respectively. We first define the operator 

Xab = (1a ® Pb) ~ Pab (18) 
It is easy to check that Xab is positive semi-definite if pab is separable. Indeed, in such a case we have 

x AB =j2 w i(( i A-p { i ) )®p i B ) )>° (19) 




since a sum of positive operators is a positive operator. Now, we can use the fact that, if X and Y are two Hcrmitian 
operators such that X > Y > then logX > logF, as implied by Lowner's theorem pj^ . (Note that the converse 
is not true.) As a consequence, using X = 1 A <8> Pb and Y = pab, we conclude that Xab > implies gab > 0. □ 

Corollary 1: Any separable bipartite state is such that pa\b < !• 

Since we have pa\b — ex p( — gab), Theorem 2 shows indeed that no eigenvalue of the conditional amplitude operator 
exceeds 1 for a separable state, as was conjectured in Ref. I^-Q]. This yields a simple necessary (but not sufficient) 
condition for separability. The classical analog of this property is that — \ogp(a\b) > 0, Va, b. The latter inequality 
simply results from the fact that p(a\b) = p(a, b)/p{b) < 1, Va, 6, as p(b) = ^2 a p(a, b) and p(a, b) > 0. 

Corollary 2: The conditional von Neumann entropy S(A\B) is non-negative for a separable bipartite state. 

Since we have S(A\B) = Tr[p AB gab], this simply follows from the fact that Tr[AF] > if X,Y > 0. Thus, 
the non-negativity of the conditional entropy is another (weaker) necessary condition for separability [^j-Q. This 
has been shown in general for "a-entropies" in Ref. It can also be related to the non- violation of entropic Bell 
inequalities ]l2[ ). 

Note that Corollary 2 can also be obtained by using the concavity of S(A\B) in a convex combination of pab's, 

S(A\B)=S(pab)-S(p b )>J2^(s(p { 2b)~S(p%)) if p AB = (20) 

i i 

a property related to the strong subadditivity of quantum entropies [^3|. Using the fact that, for a separable state, 
each term i gives S(A\B) = S(A) because p^ B — p^ ® p B (i.e., A and B are independent), we obtain 

S(A\B)>J2 w i S (pA)>° ( 21 ) 

i 

Note that a negative conditional von Neumann entropy S(A\B) necessarily implies that an eigenvalue of pa\b exceeds 
1, but the converse is not true. Thus, weak inseparability (in the sense that S(A\B) > despite the inseparability of 
Pab) may be revealed by the spectrum of Pa\b- 

Theorem 3: There exist inseparable bipartite states pab such that the operator gab is positive semi-definite; 
consequently, gab > (or pa\b 1) is not a sufficient condition for separability. 

Let us consider a bipartite system AB characterized by pab-, which we extend with another system A'B' in the 
state pa'b 1 - The joint system is then characterized by a density operator of the product form 

PAA'-BB' = PAB ® PA'B' (22) 

We first calculate the conditional amplitude operator of the joint system [AA' conditional on BB') 

PAA'\BB> = evp[\0g PAA'-BB' - log(l AA > ® PBB')} (23) 
where the reduced density operator describing BB 1 is 

PBB' = Tr^A' [PAA'-BB'} = PB ® PB' (24) 



J Here and below, the notation X > Y means that X — Y is a positive semi-definite operator. 



5 



Using the identity log(X (8 Y) = log A" ® 1 + 1 (g) logF for operators X, Y > as well as its exponential, i.e., 
exp X ® exp y = exp(A" eg) 1 + 1 ® y), we obtain 

PAA>\BB> = exp[log PAB <8 lA'B' + lAB <8>logpA'B' - 1a ® log Pb ® Ia'B' - lAB ® 1a' <8> log/OS'] 

= exp[(log p AB - 1a <8> logps) ® Ia'B' + Iab <8> (log pa'B' - 1a' ® logos')] 

= exp[logpAB - 1a ® logps] ® exp [logos' - 1a' ® logos'] (25) 

Thus, we have 

PAA'IBB' = PA|B ® PA'|B' (26) 

which parallels the classical relation p(aa'\bb') — p(a\b)p(a'\b') if AB and A'B' are independent bipartite systems, 
that is, if p(a, a'; b, b') = p(a, b)p{a! , 6'). In particular, we have 

0"(PAA'|BB') = cr(pA|s) ® O'(pA'IB') (27) 

where cr(jo) stands for the spectrum of p. 

Now, let us assume that AB is an inseparable system with gab ^ or pa\b ^ !• I n other words, the operator Pa|b 
admits an eigenvalue that exceeds 1. Assume also that A'B' corresponds to two independent systems in a product 
state, that is, pa'B' — PA' ® Pb>- The resulting conditional amplitude operator for A'B' is then pa'\b' — PA' <8> 1b', 
just like its classical counterpart p(a\b) — p(a) if p(a, b) = p(a)p(b). Obviously, we then have pa'\b> < 1> a s expected 
since A'B' is separable. According to Eq. (|27|), the eigenvalues of /9aa'|bb' are the pairwise products of eigenvalues 
°f Pa\b with eigenvalues of Pa'\b' ■ Therefore, it is easy to find a system A'B' with eigenvalues of pa'\b' small enough 
so that the product of any of them with an unclassical (> 1) eigenvalue of pa\b results in eigenvalues of paa>\bb> 
that are all < 1. The extended system is then characterized by (Jaa'-.bb 1 > or paa'\bb' < 1, while it obviously 
contains an inseparable component AB. Such a dilution of inseparability (or entanglement) is always achievable with 
a system A'B' that is large enough and maximally disordered (i.e., pa'B' ~ 1a' <8> 1b')- Consequently, the condition 
that <J ab > or p A \s < 1 cannot be a sufficient condition for separability. □ 

Remark 1: Eq. ( p7| ) implies that, if AS and j4'S' are inseparable systems with pa\b ^ 1 arL d pa'\b' ^ 1> then the 
inseparability of the joint system is revealed by Paa'\bb' ^ 1- 

Remark 2: While Xab > is a sufficient separability condition for 2x2 and 2x3 systems (cf. Section III), it cannot 
be concluded that cab > or pa\b < 1 is also sufficient in these cases, as the converse of Lowner's theorem does not 
hold. Interestingly, numerical evidence reveals that only very few inseparable states of two qubits with pa\b 1 can 
be found. 

III. NECESSARY CONDITION FOR SEPARABILITY OF MIXED STATES 

A. Bipartite system of arbitrary dimension 

As we saw in the previous section, Theorem 2 results in a simple necessary condition for the separability of mixed 
states based on Eq. , which does not require the calculation of the conditional amplitude operator (although it is 
related to it). 

Definition 3: Define a linear map A which maps Hcrmitian operators on Hab into Hermitian operators on TLab'- 



A : pab * ^ab — 1a ® Pb ~ Pab with p B = Tr A [pab] (28) 

It commutes with a unitary transformation acting independently on A and B. Indeed, according to Eq. (|l5|), if pab 
undergoes a unitary transformation of the product form JJ a ® Ub, then 

\ab^\'ab = (Ua®U b )\ab(U a ®U b ) (29) 



i.e., Xab transforms just like pab- Therefore, the spectrum of Xab is invariant under a Ua ®Ub isomorphism on 
Pab, as expected. 

Theorem 4: A necessary condition for the separability of the state pab of a bipartite system AB is that it is mapped 
by A into a non-negative operator, i.e., ApAB > 0. 
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The proof of this theorem is contained in the proof of Theorem 2 where we used the fact that the map A reveals 
non-separability: if Xab 0, then pab is inseparable. Conversely, any separable state pab is mapped into a positive 
semi-definite operator Xab- Moreover, it is easy to see that the map A conserves separability since it is linear and 
maps product states into product states: if pab is separable, then Xab > is also separable. Let us now calculate 
the partial traces of the operator Xab- 

X A = Tt b [Xab] = 1 A - pa (30) 
X B = Tt a [X A b] = (<1a - (31) 

where <1a is the dimension of TLa- This shows that the map A does not preserve the trace in general. Indeed, the 
total trace is scaled by an integer factor under A, that is, Tr[A,4s] = (cIa — 1)Tt[pab]- Thus, A is trace-preserving only 
in the special case where A is a 2-state system (i.e., d A = 2). It is also interesting to note that the map A is always 
reversible, the inverse map being given by 

A" 1 : Xab -> (d A - 1) _1 (1a <8> X B ) - X AB = Pab (32) 

where Xb is defined as above. Notice that the map A is equal to its inverse A -1 only in the case where cLa = 2. This 
implies that, if Xab is separable, then A -1 : Xab —> Pab > 0. (The fact that the inverse map reveals inseparability 
is true in this case only.) 

The separability condition based on Theorem 4 is illustrated in Appendix A, where we consider several separable 
and inseparable states. It appears that Xab > results in the same condition as Peres' in the case of two quantum 
bits, in which case it is sufficient (see Theorem 9). For larger dimensions, it is only necessary. 

Remark 1: More generally, following the approach of Horodecki et al. ||, the map A can be written as the tensor 
product of a positive linear map T and the identity, that is 

A = T ® / with r : p -> (Trp) - p (33) 

where T acts on Hermitian operators on Ha and the identity acts on operators on H.b- Since T is a positive map 
(i.e., it maps positive operators into positive operators), A = T <g> / maps separable states into positive operators 
It therefore results in a necessary condition for separability, according to Theorem 4. The map T commutes with an 
arbitrary unitary transformation U, that is 

r(C/p[/ t ) = U(T P )Ul (34) 

which makes the separability condition based onA = r® / independent on the basis chosen for A and B. In the same 
manner, the inverse map A -1 can be written as 

K- l =T- l ®I with : p -> — - p (35) 

where d is the dimension of the Hilbert space of p. Note that T _1 is not a positive map for d > 2, so that A -1 is in 
general useless as far as detecting inseparability is concerned. This emphasizes that the separability criterion based 
on Theorem 4 is quite special in 2 dimensions (e.g., for a spin-1/2 particle or a quantum bit), as will be studied in 
detail later on. Let us only mention here that the map T applied to a two-dimensional system amounts to a spin flip 
(this can be interpreted as time reversal). In order to see this, let us write an arbitrary state of a two-dimensional 
system as 

P=\{l+r-a) (36) 

where a represent the three Pauli matrices and r = Tr (pa) is a real vector in the Bloch sphere (of radius 1). The 
vector r describes the statistics of measurements on the system, as, for example, the quantum expectation value of 
the spin component along an axis defined by the vector v is 

Tr p(v ■ (?) = (v, r) (37) 

Using Eq. (pq), it is straightforward to check that 
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Tp = l-p=~(l-?-cr) (38) 

Thus, the map Y performs a spin-flip, or, equivalently, a parity transformation r — ► — r on the vector characterizing 
the system in the Bloch sphere. 

Remark 2: It is interesting to consider the classical analog of the maps Y and A = T ® J to gain some insight into 
their physical meaning. First, the map T applied to a classical probability distribution pi (diagonal p) corresponds to 
the transformation: 

Pi -> Qi = ^Pk -Pi (39) 

k 

(Obviously, qj > is not normalized except for a binary distribution.) The classical analog of the map A = Y ® I is 
thus 



Pij -> qij 



= ~ Pi\i P'-i = Pi ~ Pi i ( 40 ) 



Since puj is a probability distribution in i, we always have 1 — puj > so that qij > and the separability criterion 
is fulfilled. This emphasizes that quantum inseparability ("qij < 0") may be viewed as resulting from a conditional 
probability that exceeds 1 (more precisely, an eigenvalue of pa\b which exceeds 1). 

Remark 3: It is instructive to consider the action of the map A = Y ® I on a product state \ip) = \a) <Z> \b). Using 
Pab — P a ® Pb with P a = l a )( a l an d Pb — \b)(b\, it is easy to see that 

Xab = Pa ® Pb (41) 

where P^ = Y(\a)(a\) = 1 A - \a)(a\ is the projector on the subspace orthogonal to \a). In the special case where A is 
a 2-state system (d,A = 2), P^ is a rank-one projector as the total trace is preserved. We have then P^ = |a ± )(a ± |, 
where |a ) is a state vector orthogonal to \a) obtained by applying a complex conjugation on the components of \a) 
followed by a rotation by an angle 7r about the y-axis. (Note that it can be shown that it is impossible to construct 
a state la 2 -) that is orthogonal to an arbitrary state |o) by applying a unitary transformation alone.) Indeed, an 
arbitrary state 

\a) = a\0) + (42) 

with \a\ 2 + \(3\ 2 = 1, is transformed into 

|a ± }--/3*|0} + a *|l) (43) 
by applying to the state vector a*|0) + /3*\1) the rotation 

U y = exp(-iTray/2) = -i(7 y = a x cr z = (i q 1 ( 44 ) 

that is, a bit- and phase-flip. The transformed state {a- 1 ) is such that (a |a) = and |a^ L )(a J -| — 1a — \a)(a\, 
as expected. Consequently, the Y map (i.e., the complex conjugation followed by U y rotation) is an antiunitary 
operation on state vectors in the 2-dimensional Hubert space 'H.a- Indeed, for any two state vectors \a) and \a), we 
have (a -1 " la -1 ") = (a\a)* . Since time-reversal T is known to be an antiunitary operator given by complex conjugation K 
followed by a rotation 1Z V of angle ir about the y-axis, i.e., T = !Z y K (see, e.g., |T^|), we conclude that the map Y is 
time-reversal when d,A — 2. [This is also clear from Eq. (^8|).] Consequently, when d,A = 2, A = Y <S> I corresponds to 
physical time- reversal on the subsystem A (while leaving the subsystem B unchanged). Such a link between "local" 
time-reversal T ® I and separability has recently been pointed out by Sanpera et al. ||. The discussion of the case 
d-A = ds = 2 will be continued in Section Illb. 

Lemma 3: If pab is a separable state, then Xab is a separable operator obtained by replacing the states \a) in Ha 
by projectors P^ orthogonal to them. 

Let us assume that the bipartite system AB is separable, that is 
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PAB = 2Z Wi (l a *)( a *l ® IM(^I) 



(45) 



where the |a») <8> are pure product states [using the spectral decomposition of and p^ , it is easy to rewrite 
Eq. (0) into this form]. As a result of Remark 3, we see that pab it is mapped by A into a separable operator of the 
form 

\AB=^w i {Pt i ®\bi){bi\) □ (46) 

i 

The operator Xab is a unit-trace operator in the case <1a — 2 since each component pure state \a) (g> \b) is mapped 
into a pure product state, {a,- 1 ) <E> \b), in which case it simply reads 

X AB =J2 w i(\^)(^\^\ b i)(bi\) (47) 

i 

This expression results in a simple necessary condition for separability (distinct from Xab > 0), inspired from the 
condition recently proposed by Horodecki as follows. 

Theorem 5: A necessary separability condition for the bipartite state pab is that its support can be spanned by a 
set of product states which are such that the corresponding product states obtained by applying T to the state vector 
in the A space span the support of Xab = ^-Pab ■ 

We only consider this condition in the case where cLa — 2. The central point is to note that, if pab is separable 
as in Eq. (pTjf), then the ensemble of product states |a<) (g) \bi) span the entire support of pab @j- (Conversely, any 
state \a,i) ® \bi) must belong to the support of pab and cannot have a non- vanishing component orthogonal to it.) 
Using Lemma 3, we see that the ensemble of states \ aj~) ® \bi) span the entire support of the corresponding separable 



state Xab obtained by applying A on pab [cf. Eq. (47)]. (Also, any state \af) ® \bi) cannot be outside the support 
of Xab-) This results in another necessary condition for separability which can be stated as follows. If a state pab 
is separable, then it must be possible to span its support by a set of product states \a)\b) which are such that their 
image (i.e., the product states obtained by rotating the complex conjugate of state vector \a) in the A space by an 
angle 7r about the y-axis while leaving the state vector \b) in the B space unchanged) span the support of the mapped 
state Xab — ^-Pab- n 

Definition 4: Two additional maps from operators on TLab to operators on TLab can be defined: the dual map 

A : pab -> Xab = Pa ® Is - Pab (48) 

and the symmetric map 

M : pab Pab — 1a ® 1_b — Pa ® Is — l/i ® Pb + Pab (49) 

where p A = Tr b[pab] and p B = Tr a[pab] 

The map A which we considered until now is related to the conditional amplitude operator of A conditionally on 
B, that is pa\b- Of course, a similar linear map can be defined using the amplitude operator pb\a, and exactly the 



same conclusions follow. This is the dual map A defined in Eq. (48). It is trace- preserving and self- inverse in the 
special case where c?b = 2. It can obviously be written as the tensor product A = I ® T, where the map T now acts 
on operators on Hb, and therefore commutes with a Ua ®Ub isomorphism. Since T is positive, A maps separable 
states into positive (separable) operators, which results in another separability condition, i.e., Xab > 0. As we will 
see in Section Illb, the operators Xab and Xab can be shown to have the same spectrum when dA — ds = 2 (i.e., for 
two 2-state systems), in which case they result in the same separability condition. However, this property does not 
hold in larger dimensions, i.e., Xab and Xab do not have the same spectrum in general (see Appendix A). 

We can also construct another linear map by cascading A and A (the order is irrelevant), which results in the 



symmetric map M = A A = T ® T defined in Eq. (49). Any separable pab is mapped by M into a separable operator 
Pab > 0, as expected. The symmetric map also commutes with a Ua ®Ub isomorphism, 

M ((U A ® U b )pab{u\ ® C/jj)) - (U A ® U B ){Mp AB ){U A ® U%) (50) 

so that the spectrum of pab — M pab is invariant under local transformations on pab- It is also reversible, its inverse 
map M _1 = r _1 ® r _1 being given by 
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M 1 : [lab -> 1a ® 1b - (d B - 1) 1 {pa ® Is) - (d A - 1) x (1a ® Mb) + = Pab (51) 

where pa = Tr s [pab] = (ds - a~ Pa), Pb = T^a\mab] = (d A - — Pb)j and cU (d B ) is the dimension of Ha 

{TLb)- As expected, this map is trace-preserving and self-inverse in the special case where d A — d B — 2. It corresponds 
then to a time-reversal operation applied to the joint system AB. Note that, in this case, the symmetric map M by 
itself is not useful as far as revealing inseparability is concerned since it is positive, i.e., Mp AB > 0. Therefore, all 
inseparable states of two quantum bits are mapped into positive operators pab just as separable states. Still, M is 
important when analyzing the separability of two quantum bits as it is also equivalent to the conjugation operation 
in the "magic" basis introduced by Hill and Wootters |l(]] (see below, Theorem 10). The case of two quantum bits 
will be studied later on. Whether the positivity of M holds in arbitrary dimensions is not known. 

Theorem 6: The criterion that A maps the state pab into a positive semi-definite operator \ A b is not a sufficient 
condition for the separability of pab- 

The fact that Aab > is necessary for separability (Theorem 4) was proven before. In order to prove that it is not 
sufficient, we will show that it is possible to have an inseparable system with Aab > 0, i.e., such that its inseparability 
is not revealed by the map A. As before, we build an inseparable system p AB by extending an inseparable component 
with a separable one. Let us consider an inseparable system A'B' with \ A i B i 2 0. We extend A'B' with a separable 
system A" B" , and apply the separability criterion to the joint system AB where A = A' A" and B = B'B" . Since the 
joint system is characterized by the density operator p AB = pa'b> <& PA"B" , its associated operator under the map A 
is given by 

\ AB = ApAB = 1a ® Pb - PAB 

= (1 A > ® PB') ® (lA" &> PB") ~ PA'B' <S> PA"B" (52) 

Thus, using the operators \ a >b> = ApA'B' = 1 A > ® Pb' - PA'B' and \ A "B" = Ap A n B ,, = 1a" ® Pb» - PA»B" 
corresponding to A applied to both component systems, we obtain the simple expression 

Aab = Aa'B' ® Aa"B" + Aa'B' ® p A "B" + Pa>b> ® Aa"B" (53) 

with Aa'B' and Aa"B" > (since A'B" is separable). Thus, the sum of the first two terms on the right-hand 
side of Eq. ( |53] ) is not positive semi-definite, while only the third one is. Therefore, Eq. (|5^) does not guarantee that 
Aab > even though the composite system AB contains an inseparable component as Aa'B' ^ 0. □ 

Note that, if both components are inseparable systems that have Aa'B' ~£. and \a"B" It 0) then Aab ^ is not 
necessarily true, so that the inseparability of the joint system AB is not always revealed by the map A. [This property 
contrasts with Eq. (p6f).] Conversely, Eq. (|53| ) implies that, if both components have Aa'B' > and Aa"B" > 0, then 
Aab > 0. The example of weakly inseparable states with a positive partial transpose (see Ref. 0) is treated in 
Appendix A, to illustrate that Aab > is not a sufficient condition in general. 

Remark: The mechanism of dilution of inseparability can be understood by examining the action of the map T on 
product states. Indeed, when applying the map A = T ® / on the state pab — PA'B' ® PA"B"i T acts on the state 
PA' ® PA" {B and B' are left unchanged by /). Let us consider a density operator of the product form p = p' ® p" . 
Since we have Tr(p) = Tr(p')Tr(p"), we see that it is mapped to 

T(p' ® p") = Tr(p')Tr(p") - p' ® p" 

= [Tr(p') - p 1 ] ® [Tr(p") - /] + Tr(p') ® p" + p' ® Tr(p") - 2p' g> p" 

= Tpf ® Tp" + Tp' ® p" +p'(g> Tp" (54) 

which implies the relation 

r = r' ® r" + r' ® /" + /' ® r" (55) 

where V (or V") stands for the same map but acting on the subspace of p' (or p") while T acts on the joint space. 
Using the same notation for the map A (i.e., A' acts on the subspace of A'B' while A" acts on the subspace of A"B"), 
the latter equation gives 

A = T $ I = A' <gi A" + A' <g> I" + I' ® A" (56) 

which results in Eq. (|53|) . The same reasoning can be applied to the dual map A = / ® T and to the symmetric map 
M = r ® T. Thus, even if the maps A' and A" reveal inseparability by themselves, the combined map, Eq. (p6|), is 
not guaranteed to do so because the non-positivity of (A' ® A")p can be masked by the two following terms. 
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B. Special case of two 2-dimensional systems 



Lemma 4: The map T acting on the state of a two-dimensional system corresponds to time-reversal, and is therefore 
equivalent to applying the complex conjugation operator K followed by a rotation !Z y by an angle tt about the y-axis, 
that is, T = Tl y K. 

Let us start by summarizing the action of the map r : p — ► (Trp) — p on the density operator p characterizing a 
2-dimensional system (e.g., a quantum bit). Since p can be written as a linear combination of the unit matrix and 
the three Pauli matrices a with real coefficients, it is sufficient to consider the action of T on these (Hermitian) basis 
matrices. It is straightforward to check that T is an antiunitary operator that leaves the unit matrix unchanged and 
flips the sign of the Pauli matrices <J x .y,z, 



(57) 



so that, in the Bloch-sphere picture, we obtain 



p = i(l + r*. a) -> Tp = 1 - p = ~(1 -r-a) (58) 

where r is the Bloch vector. Thus, the map T flips a (or performs a parity transformation on the Bloch vector r — ■> — f), 
and corresponds to time-reversal. It can therefore be decomposed into complex conjugation K followed by a rotation 
lZ y of an angle n about the y-axis, that is F = T = lZ y K Jlifj . □ 

Remark: The complex conjugation operator K (or equivalcntly the transposition, as we deal with Hermitian opera- 
tors) corresponds to an antiunitary operator which acts on the four basis matrices as: 

1 > 1 <7 X > G x (T y ► -Oy G z — ► <T Z (59) 

(Remember that it is enough to consider the action of K on the basis matrices as the coefficients are real.) Also, lZ y 
is a unitary operation characterized by the unitary matrix U y = exp(— may/2) — ~io y = a x a z which maps p into 
UypU^ — (Typcfy, so that the basis matrices are transformed according to 

i ^-y ^-v ^-y inc\\ 

— ► 1 a x — > ~g x (T y — > <j y a z — ► -<r z (60) 



It is straightforward to check, using Eqs. (|57|), ( p9|) and (|6C|), that T is the product of K and lZ y . (It is a general 
property of an antiunitary transformation that it can be written as the product of a unitary transformation and a 
fixed antiunitary operator such as time-reversal.) 

In short, one can see that, if the system is in a (pure or mixed) state given by Eq. (36), then 

u vP* U v = VyP* a y = \( l + a v(r- o*)<J y ) = -(1 - f ■ a) = Tp (61) 

where we have used the fact that r is a real vector and that a y aa y = — a*. This generalizes what was shown earlier 
for pure states, namely that if \a) = a|0) + /3|1) and \a^) = U y (a \0) + = — j9*|0) + a*|l), then we have 

\a ± ){a x \=T(\a)(a\) (62) 

Corollary 1: When writing the Hilbert-Schmidt decomposition of pab for two 2-state systems, the map A = T ® I 
corresponds to a sign-flip of the Pauli matrices acting on A while leaving the sign of those acting on B unchanged. 

Let us consider the Hilbert-Schmidt decomposition of an arbitrary state of AB (e.g., two quantum bits or spin-1/2 
particles) M: 



Pab — j I 1a ® 1b + r ■ oa ® Is + 1a <E> s ■ a B + ^ t n ,m c r i™ ) ® cr^ n) j (63) 



\ m,n— 1 / 

where a A "^ and stand for the Pauli matrices (with n = 1,2, 3) in the A and B space, respectively. Eq. ( |6^ ) 
depends on 15 real parameters, the two 3-dimensional vectors r and s, and the 3x3 real matrix t n ^ m . The vectors r 
and s correspond to the state of A and B in the Bloch sphere since we have 
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PA = Tt b [pab] = + f ■ a a) 



Pb = Ti'a[pab] = 2 (is 



s ■ o B ) 



(64) 
(65) 



They characterize the reduced systems A and B, that is the local (marginal) statistics of any observable on A or B. 
The matrix t niTn = Tr[pA B (o~A Kxrji^)] describes the joint statistics of A and B, and yields for example a very simple 
expression for the correlation between the measured spin components along two axes (defined by the vectors a and 

by. 



Tr 



(3, tb) 



p(a- &A®b- as) 

It is checked by straightforward calculation that the map A = T ® I simply flips the sign of the terms in 3a'- 
Aab = 1a <8 pB - PAB 



(66) 




r ■ a a 



s ■ a B 



» 



» 



(67) 



Corollary 2: The dual map A = I ® T flips the sign of the Pauli matrices acting on B while leaving the sign of those 
acting on A unchanged. The action of the symmetric map M on the Hilbert-Schmidt decomposition of pab is to flip 
the sign of the Pauli matrices acting on both A and B: 

Pab = 1a ® Is - Pa ® Is - 1a ® + pab 

= - I 1^ ® 1 B - r- CTA ® Is - lA ® S ■ G B + tn,m ® cr^ 1 ^ J (68) 

\ m,n— 1 / 

This operation corresponds to time-reversal applied to ^4 and B simultaneously, and can be shown to be equivalent 
to the conjugation operation in the "magic" basis introduced by Hill and Wootters Jl0| (see below). 

It is worth noting that the set of states that remain invariant under the symmetric map M are mixtures of generalized 
Bell states, the latter being defined as the states obtained by applying any local transformation to the four Bell states. 
These states are also called "T-states" by Horodecki et al. ||, and are such that the entropy of A and B is maximal, 
that is S(pa) = S(p B ) = 1. (The only pure states in this set are the fully entangled states of two qubits, i.e., the 
generalized Bell states.) Thus, in particular, the (generalized) Bell states are left unchanged by the action of M. In 
contrast, a (separable) product state pa® Pb is mapped into the distinct state p-ab = (1a — Pa) ® (1_b — Pb)- Because 
of this property, pab by itself is uninteresting as far as revealing inseparability is concerned, as mentioned earlier. 

Theorem 7: For two-dimensional systems A and B, the map M conserves the spectrum, so that the separability 
criteria resulting from the map A and its dual A are equivalent. 

As we have seen before, the map T acting on subsystem A (or B) amounts to performing a complex conjugation 
operation K followed by a rotation !Z y defined by U y = exp(— iTtt7 y /2) = —iu y (i.e., a bit- and phase-flip). Thus, the 
symmetric map M = T <g> T is equivalent to a complex conjugation K (or transposition) of the joint density operator 
in the Hilbert space Hab, followed by a tensor product of rotations U y ®U v — ~o y ® a y . Note that, as we are dealing 
with Hermitian (density) operators, their spectrum is unchanged by K . The same is true for the rotation U y (g> U y . 
As a consequence, the operator pab — MpAB has the same spectrum as pab when oIa = d B = 2. As T is self-inverse 
(r 2 = I) when cLa = d B = 2, we have the relation I ® T = (T ® I)(T ® T) or in short A = AM. This implies that 

Apab = A [{U v ® U y )p AB {Ul ® Ul)} (69) 

which in turn results in 

Xab = (U y (8 U v )X* AB (Ul ® Ul) (70) 

as A commutes with U v ® U y and complex conjugation. Since \ab is Hermitian (just as pab), the latter expression 
shows that the spectrum of A as and Xab are identical, so that the resulting criteria for separability are equivalent. 
□ 
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Theorem 8: For two-dimensional systems A and B which have maximal reduced entropy, i.e., S(pa) = S(pb) = 1, 
the positivity of the operator ApAB results in a necessary and sufficient condition for the separability of pab ■ 

The states with f = s = Q are such that the reduced density operators are given by pa = Pb = 1/2, so that 
the reduced entropies are S(pa) = S(pb) = 1- These "T-states" are thus completely characterized by the matrix 
tn,m- It has been shown in Ref. Q that any state belonging to this set of T-states can be transformed by a unitary 
transformation of the product form Ua®Ub into a state for which t n ^ m is diagonal. As far as separability is concerned, 
we can thus restrict ourselves to the class of all states with diagonal t, since these are representative of the entire set 
of T-states (up to an Ua ® Vb isomorphism). 

The class of states with diagonal i is a convex subset of the set of T-states, and any state belonging to this 
subset can be characterized by the real vector t = (in, £22, £33) made out of the diagonal elements of t. It has been 
shown in Ref. Q that an operator pab of the form given by Eq. (|6^) with r = s = and diagonal t corresponds 
to a state (i.e., a positive unit-trace operator) if and only if the vector t belongs to a tetrahedron with vertices 
t\ = ( — 1, 1, 1), t2 = (1, —1, 1), £3 = (1, 1, —1), and ti = ( — 1, —1, —1). In other words, any state of this class can be 
represented by a point inside this tetrahedron. In this representation, the four Bell states ( I >± } = 2~ 1 / 2 (|00) ± |11)) 
and |^ ,± ) = 2~ 1 / 2 (|01) ± 1 10) ) correspond to the vertices of the tetrahedron, that is 

h : |$-><$-| 
t 2 : |* + )(* + | 
f 3 : |*+)(*+ 1 

r 4 : i*->(*- 1 

In Ref. Q , it is also shown that a state pab of this T-diagonal class is separable if and only if the vector t characterizing 
Pab belongs to an octahedron with vertices o x = (±1,0,0), o 2 = (0, ±1,0), and o^ 1 — (0,0, ±1). This results in a 
necessary and sufficient condition for separability within the class of T-states. 

Let us consider the action of A in this representation. As shown earlier, A flips the "spin" ua- Within the set 
of T-states, this amounts to changing the sign of the t n ^ m matrix, that is, to flipping the sign of the vector t for 
T-diagonal states. Therefore, the criterion for separability Xab — ^-Pab > translates, in this representation, to the 
condition that the "parity" operation on the vector t characterizing a separable state results in a positive operator 
(i.e., a legitimate state). Thus, —t must belong to the tetrahedron. It is easy to see that the set of points of the 
tetrahedron which are such that their image under parity still belongs to the tetrahedron corresponds exactly to 
the octahedron defined above. Therefore, no inseparable state exists that satisfies ApAB > 0, so that A provides a 
necessary and sufficient condition for separability within the class of T-states. □ 

Theorem 9: A bipartite system of two-dimensional components A and B characterized by an arbitrary joint density 
operator pab is separable if and only if the operator Xab — ApAB is positive semi-definite. 

It is enough to show that A is equivalent to a partial transposition up to some completely positive map (in fact, 
a unitary transformation). As already mentioned, Peres' separability criterion is based on the partial transposition 
operation, that is, on the map T ® I, where T is the standard transposition of operators on the Ha subspace. Since 
we are dealing with Hermitian operators, T® / is equivalent to the "partial conjugation" operation K ® I, where K is 
the complex conjugation operator acting on states on TCa- Note that, although K is well-defined, partial conjugation 
K <g> I is only defined for product state vectors in Hab @j- We can now use Lemma 4, i.e., T = lZ y K, together with 
the fact that any positive map II acting on density operators in a two-dimensional Hilbert space can be written as H 

n = p + n^ p T (72) 

where Ilf p and IlS? p are completely positive maps (which therefore do not reveal inseparability) . With the identifica- 
tion Ilf p = and 1T2 p = TZ y , we see that the map T can be used rather than the transposition operator T (or K ) in 
order to test the positivity of the operator resulting from applying any element of the set of maps II <S> I on pab ■ This 
is simply due to the fact that the complex conjugation operator K is unitarily equivalent to T. Thus, the reasoning 
used in Ref. Q holds here, so that KpAB > results in a necessary and sufficient condition for the separability of 

PAB- □ 

Remark 1: Since the spectrum of an operator is conserved by a unitary transformation (!Z y ), it is clear that 
the spectrum of the matrix obtained by partial transposition in subspace A, Pab, is the same as the spectrum of 
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Xab = A-pab- Therefore, testing Peres' condition or the positivity of Xab is operationally equivalent, and these 
conditions can be used interchangeably in the case of two quantum bits, as illustrated in Appendix A. Since T 
corresponds to time- reversal (on the subsystem A) , the map A = T ® / amounts to changing the arrow of time for 
subsystem A with respect to subsystem B. Such a relation between time-reversal and Peres' map has been pointed 
out previously by Sanpera et al. P] , where it was shown that the partial transposition operator is unitarily equivalent 
to "local" time-reversal. 

Remark 2: As we know that T applied to a two-dimensional system is unitarily equivalent to the transposition 
operator T, the separability condition based on A = T ® I is in fact equivalent to the one based on the partial 
transpose (T ®I)pab = Pab whenever subsystem A is two-dimensional. More precisely, Xab and p^ B have the same 
spectrum for 2 x n systems, so that the conditions arc equivalent if T is applied on the two-dimensional subsystem. 
As a consequence, the separability condition based on A is also necessary and sufficient for 2 x 3 systems, while it 
is only necessary for 2 x n systems with larger n, just as Peres' condition ||. Numerical evidence suggests that, for 
systems with d,A,dB > 2, the condition based on A (or A) is weaker than (or equivalent to) the one based on partial 
transposition. 

Theorem 10: The symmetric map M = T <g> T applied to a bipartite system of two-dimensional components (i.e., 
global time-reversal) is equivalent to complex conjugation in the "magic" basis introduced in Ref. ]l6[| . 

Since T — 1Z V K, where K denotes the conjugation operator and lZ y is a rotation characterized by U y = 
exp(— iira y /2) = —i<jy, the symmetric map M applied to the state pab of a bipartite systems results in 

Mpab = (U v ® U y )p* AB (Ul ® Ul) (73) 

where U y ®> U y = —<7 y (8> a y . Since M is antiunitary and self- inverse (Af 2 = I), it is a conjugation [ pi . It can be 
written as the complex conjugation operator if expressed in a specific basis. Let us assume that V is the unitary 
operator (in the joint space) that transforms the product states into the states { | e^) } that form this specific basis, 
that is 



|ei) = V|00> |e 2 ) = y|01) |e 3 > = V|10) |e 4 ) = V\ll) 



(74) 



We would like to show that M is equivalent to rotating the states \ei) into the product states, taking the complex 
conjugation of the density matrix (in the product basis), and then rotating the product states back to the |ej)'s: 



M Pab = V{V^p AB V)*V^ = (VV T )p* AB (VV T ^ 
where V T is the transpose of the unitary matrix V. Identifying Eqs. ([73]) and (]7j]), we obtain 



(75) 



VV 1 =Uy®Uy = -<Ty (g> O y 



/ 1 \ 
0-10 
0-100 

V i o oo/ 



(76) 



It is easy to prove that, if V is unitary, then VV T is unitary and symmetric (but not necessarily Hermitian). In order 
to find a solution for V that satisfies Eq. (|76|), we first diagonalize the matrix o~ y ® a y . Consider the unitary matrix 



(77) 



W = exp ( -^j-(l - °~x) ® (1 - Ox) ) = (1 ® 1 + 1 ® o x + a x ® 1 - a x ® cr x )/2 



This matrix is self-inverse, that is, W 2 = 1, so that it is Hermitian (W^ = W). It can easily be shown that it 
diagonalizes o y ® o~ yi i.e., 



W(a y (gi a y )W = a z ® a z 

Replacing V in Eq. (|76|) by the product of a diagonal matrix D and W , that is V — WD, we obtain 



(78) 



DD 1 



-W(0-y®(Jy)W 



-a z ®a z = 



I -1 \ 

10 

10 

\ -1 / 



(79) 
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which implies that 
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This yields a (non-unique) solution for the unitary matrix V = WD which defines the basis {|ej}}. It is worth noticing 
that the matrix 



W 



( 1 1 1 -\\ 
11-11 
1-111 

V-i i i i / 



(81) 



transforms the product states into the four maximally entangled states which are obtained by applying a local 
transformation H <£> 1 on the four Bell states, i.e., 

W\QQ) = (|00) + |01) + |10) - |ll))/2 = (H ® 1)|#+) 
W|01) = (|00) + |01) - |10) + |ll))/2 = (H ® 1)|*+) 
W\1Q) = (|00) - |01) + |10) + |ll})/2 = (if ® 1)|$-) 

W\ll) = (—|00) + |01> + |10) + |ll»/2 = (H® l)|*-> (82) 

where H is the Hadamard transform. (As a matter of fact, the unitary transformation W corresponds simply to 
a controlled-NOT gate where the control is in the dual basis {|0) + |1), |0) — |1)} rather than the standard basis.) 
Therefore, the unitary transformation V = WD is such that the product states are rotated into the four generalized 
Bell states with the appropriate phases 

|ei) = V|00) = ±i(H® 1)|$+) 
\e 2 ) = V\0l) = ±1{H (g> 
|e3> = V|10) = ±l(fl"(8l)|*-) 
|e 4 ) = V|ll) = ±i(H®l)\y-) 

(83) 

These states |e») are therefore equivalent, up to a local change of basis H ® 1 and a phase i that are irrelevant here, 
to the "magic" states introduced in Ref. Jl^]. (Any four states obtained from the |ej)'s up to an overall phase and a 
unitary transformation acting locally on each quantum bit are legitimate "magic" states.) This implies that, when 
expressed in this basis, the symmetric map M = T ® T reduces the the complex conjugation operation that was used 
in the context of the calculation of the entropy of formation of a pair of quantum bits (see Refs. 00). □ 



IV. CONCLUSION 



Given a bipartite system characterized by a density operator pab , one can define a conditional amplitude operator 
Pa\b ( a positive Hermitian operator defined on the range of pab) which plays the role of a conditional probability 
in quantum information theory. Specifically, it can be used to define a conditional von Neumann entropy, S(A\B) = 
— Tt[pab logpA|s]j m perfect analogy with the conditional Shannon entropy. The quantum counterpart of many 
classical properties also holds: i) pa\b is defined on the support of pab, so that S(A\B) is well-defined; ii) S(A\B) = 
S(AB)-S(B); iii) p A \ B = Pa®~Lb if A and B are independent; iv) p A A>\BB> = Pa\b®PA'\B' '^Paa'-bb' = Pab®Pa'b>\ 
v ) Pa\b transforms as (U a ® Ub)pa\b{U\ ® U B ) when performing a local unitary transformation Ua <8> Ub on pab, 
so that its spectrum and therefore S(A\B) are invariant under such transformations on pab- 

The main non-classical feature that appears when dealing with a quantum bipartite system rather than a classical 
one is that pa\b may have a "non-classical" spectrum, that is, eigenvalues of Pa\b m ay exceed 1, which in turn 
implies that S(A\B) can be negative. More specifically, we have shown that pab < 1 for any separable state, which 
also straightforwardly implies S(A\B) > 0. Therefore, a necessary condition for separability is that the conditional 
amplitude operator has a "classical" spectrum, or that the conditional entropy is non-negative (the latter is a weaker 
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condition). These conditions are not sufficient, since extending an inseparable state with a separable one of large 
dimension may result in a dilution of inseparability, that is, it may give rise to a state with Pa\b < 1. In other words, 
some inseparable states exist with pa\b < 1, and certainly some with S(A\B) > (even if pa\b i !)• 

These considerations can be used to define a simpler necessary condition for separability, based on the positive 
linear map V : p — » (Trp) — p. Any separable state is mapped by the tensor product of Y (acting on A) and the 
identity I (acting on B) into a positive operator. Therefore, a simple separability criterion is based on checking the 
positivity of the operator (T ® I)pab = Pb — Pab- This condition, along with the one based on the dual map 
I <8> r, can be shown to be non-sufficient for a system of arbitrary dimension. Since the map Y commutes with any 
unitary transformation, the spectrum of the operator (Y (g> I) Pab is invariant under a local unitary transformation 
Ua <E> Ub, making this condition independent of the basis in which A and B are expressed. 

In the case of a 2-dimensional system, the map Y can be shown to be the time- reversal operator, which flips the 
sign of the spin matrices (or, consequently, reverses the Bloch vector characterizing the state of the quantum bit). It 
is therefore unitarily equivalent to the transposition operator |^|| , so that the condition based on Y / is equivalent 
to the one based on Peres' partial transposition for 2 x n systems (when applying the map Y on the 2-dimensional 
subsystem) . As a consequence, it is necessary and sufficient for 2x2 and 2x3 systems while it is only necessary for 
larger systems, just as is Peres' E). Numerical evidence suggests that, for systems where dA,ds > 2, the condition 
based on Y ® I or I (g) Y is weaker than (or equivalent to) the one based on partial transposition. 

Finally, we consider the symmetric map, defined as (r <X> Y)pAB — 1a ® Is ~ PA ®^-b~ 1a ® Pb + Pab- The 
states which are left invariant under this map are mixtures of generalized Bell states (or "T-states"), which include 
the maximally entangled pure states as well as the product of two independent (unentanglcd) random bits. It can be 
seen that the map r ® Y is also related to quantum nonlocality of two quantum bits even though it does not directly 
reveal inseparability. Indeed, Y ® Y reduces to the complex conjugation in the "magic" basis that has been used in 
the context of the calculation of the entropy of formation of a pair of quantum bits (see Refs. jl0],[l7]] ) . It therefore 
might be interesting to look for a simple relation between the map Y (related to inseparability) and the entropy of 
formation. This will be the subject of further work. 

Note: Parts of Section MB of this paper are equivalent to results contained in Ref. [[l7|], which was brought to our 
attention after completion of this work. 
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APPENDIX A: EXAMPLES 



Here we consider several examples illustrating the separability criterion Xab > 0, and compare it to Peres' crite- 
rion M. Examples 1-4 deal with states of two quantum bits, and illustrate the fact that the A-criterion is necessary 
and sufficient (the spectrum of Xab is identical to the spectrum of p TA )- Examples 5-6 illustrate that the A-condition 
is not sufficient for systems in larger dimensions (3x3 and 2x4) whose partial transpose is positive (cf. Ref. H). In 
fact, the A-condition is equivalent to Peres' condition for 2 x n systems, so that it is also necessary and sufficient for 
2x3 systems while it is only necessary for larger n. 

Example 1: Consider a Werner state jjj with parameter x (0 < x < 1), that is a mixture of a fraction x of the 
singlet state l^ - ) and a random fraction (1 — x). We shall see that Xab > is equivalent to Peres' criterion, and is 
therefore sufficient in this case. Indeed, the joint density operator 



Pab =sc|* )(*" 



(1-x) 



(1(8)1) 
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1-x 
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1+x 
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1+x 
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1-x 
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is mapped by A into the operator 



(Al) 
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Xai 



/i±* 
^ 



V 



\ 




1+x 
4 / 



(A2) 



which admits three eigenvalues equal to (1 + x)/4 and a fourth equal to (1 — 3x)/4. The latter becomes negative if 
x > 1/3, so that Xab is positive semi-definite only if x < 1/3, which has been proven to be the exact threshold for 
separability (any Werner state with x < 1/3 is separable as it can be written as a mixture of product states As 
expected, the spectrum of Xab is equal to the spectrum of the partial transpose of Pab-, so that the A-condition is 
sufficient to ensure separability for Werner states. 

Example 2: Consider a mixed state that is made out of a fraction x of the entangled state l^) = a|01) + 6 1 10) , and 
fractions (1 — x)/2 of the separable product states |00) and |11) (see @). The joint density matrix is of the form 



PAB = x\ip){ip\ + 



|00}(00| 



-lll)(lll = 



/ 1=2 
/ 2 









cc|a| 2 xab* 

xa*b x\b\ 2 





° \ 





2 I 



(A3) 



with a and b satisfying \a\ 2 + \b\ 2 = 1. It is mapped by A into the matrix 



Xab = 



( A*>? 
o 



V o 





l-x 
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-xa*b 






-xab* 

l-x 
2 





(A4) 



x\a\ 2 J 



The eigenvalues of Xab are:r|a| 2 , x\b\ 2 , and (1 — x±2x\ab\) /2. This implies that pab is inseparable if a; > (l + 2|a6|)~ 1 , 
exactly as predicted by Peres using the partial transpose of pab- Since we are dealing with two qubits, this is the 
exact limit between separability and inseparability 



Example 3: In the simpler case where pab is a mixture of a fraction x of the singlet state 1^ ) and a fraction (1 — x) 
of the separable product state |00), 



PA B=a;^>(^| + (l-aj)|00>(00| = 



l-x 









\ 


x/2 


-x/2 


-x/2 


x/2 





J 



(A5) 



we obtain 



Xab = 



fx/2 \ 

x/2 

x/2 l-x 

\ x/2 J 



(A6) 



2, so 



The latter matrix admits two eigenvalues equal to x/2 and two eigenvalues equal to yl — x ± \fjl — x) 2 + : 

that its determinant is equal to — (x/2) . Thus, this state is inseparable whenever x > 0, as expected. (It is separable 
only if it is the pure product state |00).) 



Example 4: Consider the class of 2-qubit inseparable states described by Horodecki et al. 
entangled states: 



where \ipx) 
matrix 



PAB = P\ipl)(i>i\ + (1 -f>) IV'aXV'a | 
a|00) + and \ip 2 ) = a\01) +£>|10), with a, b > and satisfying \a\ 2 



\b\< 



PAB = 



(pa 2 pab\ 

(l-p)a 2 (l-p)ab 

(l-p)ab (l-p)b 2 

\ pab pb 2 J 



|^], a mixture of two 
(A7) 

1. The joint density 
(A8) 
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is mapped by A to 



XaI 



/(l-p)6 2 

pb 2 (p - l)ab 

[p— l)ab pa 2 
\ -pab 



— pab \ 



[l-p)a 2 / 



(A9) 



The latter matrix admits two eigenvalues equal to (p ± ^Jp 2 + 4a 2 6 2 (l — 2p)j j 2 and two eigenvalues equal 

2, so that its determinant is equal to — a 4 6 4 (l — 2p) 2 . 



to 



1 - p ± ^(1 - p) 2 + 4a 2 b 2 (2p - 1 
fore inseparable whenever ab ^ and p ^ 1/2, in perfect agreement with Ref. || 

Example 5: Consider the 3x3 system in a weakly inseparable state introduced by Horodecki 



This state is there- 



1 



PAB 



1 + 8a 



/ 


a 











a 











a 





a 































a 































o 



















a 











a 











a 



















a 































l+a 
2 





2 

























a 





V 


a 











a 





2 





l+a 
2 



(A10) 



/ 

where a is a parameter (a ^ 0, 1). As shown in Ref. Q, the partial transpose of this state is positive, although pab 
is inseparable, which makes the inseparability of pab undetectable using Peres' criterion. It is simple to check that 
the A-mapped operator 



1 



l + 8a 



l+3a 
2 





Vl-a 2 
2 





—a 











—a 





2a 























Vl-a 2 
2 





l+3a 
2 





























l+3a 
2 





v/l-a a 
2 











—a 











2a 











—ft 











v/l- 2 
2 





l+3a 
2 





























2a 





























2n 





—a 











—a 











2a 



(All) 



V 

is positive (with a trace equal to 2), so that A cannot reveal the inseparability of pab either. Accordingly, the 
determinant of Xab is equal to 6a 7 (l — a)(5a + 3)/(l + 8a) 9 and thus positive. Note that the dual map also yields a 
positive operator Xab (of trace 2), although the eigenvalues of Xab are distinct from those of Xab, as is its determinant 
Det(AyiB) = 24a 7 (1 — a 2 ) / (l + 8a) 9 . This example emphasizes the fact that A does not result in a sufficient separability 
condition for 3 x 3 systems, just as Peres' condition ||. 

Example 6: Following Horodecki Mf, we consider a 2 x 4 system in an inseparable state 
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(A12) 



that has a positive partial transpose, where 6 is a parameter (6 ^ 0, 1). Applying A, we see that 
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(A13) 



has eigenvalues 0, 6, 26, and ^1 + 26 ± + 2b) 2 — 26(3 4- 6) j j 2 so that it is always non-negative. Note that the 

spectrum of Xab is the same as the spectrum of the partial transpose p^s ( c f- ®)j as expected. This confirms that 
the condition based on A = T®7 and Peres' separability condition are equivalent for 2xn systems (when T is applied 
to the two-dimensional system and I to the n-dimcnsional one). In this example, applying the dual map A = I ® T 
yields a positive operator which traces to 3. 
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